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HIGH  ENERGY  ION  EXPANSION  IN 
LASER.PLASMA  INTERACTIONS 


Historically,  the  high  energy  ions  in  a laser-produced  plasma  are 
defined  as  a small  group  of  ions  transporting  a significant  fraction  of 
the  absorbed  laser  energy.  Most  expansion  models  ' indicate  that  the 
energetic  ions  are  the  direct  consequence  of  the  presence  of  high 
energy  electrons.  Most  plasma  simulations  ’ use  a single  ion  species 
to  model  the  ion  expansion.  Here  we  show  that,  although 
the  hot  electron  expansion  can  account  for  the  energy  content  of  the 
fast  ions,  a multi-ion  species  description  is  usually  required  to 
reproduce  the  measured  high  energy  ion  distributions.  Ion  energy 
distributions  measured  from  CD^  targets,  where  both  predominant 
species,  C and  D , have  the  same  charge-to-mass  ratio,  are  adequately 
represented  by  a single-ion  species  expansion.  For  CH^  targets, 
however,  a two-ion  fluid  description  is  required  to  reproduce  the 
qualitative  features  of  Che  ion  expansion. 

Typical  single-shot  energy  distributions  of  high  energy  ions 
measured  from  the  Nd- laser  irradiation  (75  psec,  ~ 10^®  W/cm^  with  an 
f/2  lens)  of  a CD^  and  CH^  target  are  shown  in  Figs,  1 and  2 
respectively.  The  gross  variation  between  the  charge  collector  “ 


oscilloscope  traces  such  as  shown  in  Figs,  lb  and  2b  are  due  mainly  to 

the  different  viewing  angle  between  the  target  normal  and  the  analyzer  ^ 

axis.  For  both  cases  the  analyzer  axis  was  at  35°  with  respect  to  the 

laser  axis  but  the  target  normal  in  Figs.  1 and  2 were  55°  and  22®  lUBUnY  CClltS 
Note:  Manuscript  submitted  November  23,  1977. 
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respectively  from  the  laser  axis  toward  the  analyzer.  The  high  energy  • 

I 

electrostatic  ion  analyzer‘s  used  has  twelve  channels  with  IC^  species 

+s  + * 

and  energy  resolution.  Since  this  analyzer  cannot  resolve  C from  D 
ions,  the  ions  were  removed  by  allowing  them  to  charge  exchange  with 
a nitrogen  gas  background  (1.5  X 10'“^  Torr)  in  the  target  chamber  before 
entering  the  ion  analyzer.  The  lowered  ionization  stages,  originally 
C^,  can  then  be  differentiated  from  d"*”  ions.  The  ion  energy 
distribution  is  finally  reconstructed  from  the  summation  of  the 
recombination  products.  For  either  CH^  or  CD^  targets  no  ionization 
stages  lower  than  and  only  a small  fraction  of  appeared  under 
good  vacuum  conditions  (S  x lO”^  Torr). 

The  multipeak  structures  on  the  H and  D energy  distributions  ^ 

were  always  observed  but  were  nonreproducible  in  detail.*^ The  main 
difference  between  the  two  targets  is  in  the  relative  behavior  of  the  * 

4ia  ^ ^ 

C ion  energy  distribution  with  respect  to  the  H or  D ion 
distribution.  For  the  CD^  target,  the  ratio  of  the  number  of  to 
ions  remains  approximately  constant  with  increasing  E/Z  (energy 
divided  by  the  charge).  For  the  CH^  target,  however,  little 
correlation  was  found  between  the  and  H peaks  but  the  average 
ratio  of  H to  ions  was  increasing  with  E/Z.  In  fact,  above  50  kev/ 

Z,  more  than  half  the  fast  ion  energy  was  transported  by  ions. 

Figure  1 suggests  that  the  plasma  expansion  from  the  CD^  target 
behaves  like  a single  ion  species  during  the  acceleration  phase  of  the  ^ 

expansion.  The  relative  behavior  of  and  ions  should  therefore  , 

be  predictable  from  a single  particle  model.  From  the  equation  of 
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I 


motion  of  an  ion  accelerated  in  an  electric  field  S the  rate  of  change 
of  ion  kinetic  energy  ^ given  by 

• (« 

i ip 

where  and  are  respectively  the  atomic  number  and  charge  state 
of  the  ion.  From  Eq.  1,  only  the  ion  species  with  the  same  A/Z  acquire 
the  same  E/Z  regardless  of  the  electric  field  spatial  or  temporal 
dependence.  Both  and  D have  an  A/Z  of  2.  Figure  1 also  shows 

4b  f 

that  CD  targets  yield  C and  D ions  with  the  same  E/Z  behavior. 

2 

This  leads  to  the  conclusion  that  most  of  the  carbon  ions  were  already 
fully  ionized  before  the  acceleration  phase.  A mixture  of  ionization 
stages  lower  than  could  not  give  the  same  E/Z  for  and  D ions 
if  the  carbon  ions  were  stripped  during  or  after  acceleration.  A CD^ 
target  can  therefore  be  adequately  represented  by  a single  ion  species. 

For  a CH^  target,  Eq.  1 can  be  used  for  the  case  of  different  ion 
species  accelerated  through  a static  potential  sheet  for  different 
acceleration  times,  i.e.,  t « (Aj^/Z^^)*.  Under  this  assumption, 

4b  + 

C and  H should  still  have  the  same  E/Z  behavior.  Figure  2, 

however,  contradicts  the  static  potential  assumption  since  an  important 
+ 

fraction  of  the  H ions  is  accelerated  to  an  E/Z  larger  than  the  C 
ions.  A second  approach  is  to  assume  that  the  expansion  of  the  faster 
Ions  decreases  the  electric  field  strength  such  that  the 
acceleration  time  and  the  electric  field  are  essentially  the  same  for 
both  and  H ions.  The  final  energy  relationship  between  and 
ions  is  then  idealy  given  by 
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An  initial  ion  temperature^  of  a few  hundred  eV  also  makes  the 
viscosity  term  (ocT^^)  negligible  with  respect  to  the  ambipolar 
electric  field.®’®  The  ion  temperature,  although  high  enough  to 
neglect  viscosity,  remains  relatively  small  compared  to  the  electron 
temperature,  so  that  the  ion  pressure  can  also  be  neglected. 

An  assumption  about  the  electron  temperature  and  the  heat  flow  is 
required  to  close  the  moment  equations.  For  the  case  presented  in 
Fig.  3,  a uniform  electron  temperature  throughout  the  expansion  region 
was  assumed,  i.e.,  the  heat  is  allowed  to  flow  without  inhibition.  The 
left  boundary  in  Fig.  5a  is  an  impenetrable  wall.  The  total  energy 
is  then  conserved  by  reducing  the  electron  temperature  according  to  the 
rate  of  change  of  ion  kinetic  energy.  Other  heat  flow  assumptions  such 
as  an  adiabatic  expansion^  or  a strongly  inhibited  heat  transport 
(flux  limiter)^  have  also  been  used,  yielding  no  fundamental  differences 
in  the  qualitative  features  that  will  be  discussed  below.  Our  two- 
fluid  model  is  therefore  not  a strong  test  of  the  validity  of  the 
isothermal  expansion  assumption.®’^® 

The  set  of  fluid  equations  has  been  solved  numerically  using  an  FCT 
algorithm^^  on  a sliding-zone  grid.  Figure  5a  shows  the  evolution  of 
the  ion  density  profiles  for  a CH^  plasma  after  5.6  t (where  t is  the 
density  gradient  scale  length  divided  by  the  hydrogen  ion  sound  speed). 
As  can  be  seen  from  Fig.  5b-c,  most  of  the  ion  expansion  energy  is 
contained  in  a small  fraction  of  the  ions  with  energies  higher  than  the 
initial  electron  temperature.  About  75'*'  (Fig.  5c)  of  the  electron 
thermal  energy  remains  after  5.6  r but,  because  of  the  much  weaker 
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density  gradient,  ion  acceleration  by  the  ambipolar  electric  field  is 


greatly  reduced.  The  ion  acceleration  time  is  therefore  approximately 


the  same  for  both  C ° and  H ions.  The  remaining  thermal  energy  will 


be  dissipated  via  channels  other  than  fast  ion  production 


The  asymptotic  energy  distribution,  calculated  from  Fig.  3,  is 
shown  in  Fig.  k assuming  an  initial  electron  temperature  of  50  keV.  The 


ratio  of  H to  C ° ions  is  increasing  with  ion  energy,  in  qualitative 


agreement  with  the  measured  energy  distributions.  Other  initial  density 


profiles  have  also  been  tried  in  the  model  which  affected  the  detailed 


shape  of  the  asymptotic  ion  energy  distribution.  However,  in  all  cases 


a significant  fraction  of  the  H ions  was  always  observed  at  higher  E/Z 


than  the  ions.  From  Fig.  3,  one  can  see  that,  although  the  ion 


acceleration  time  was  about  the  same  for  both  ion  species,  the  electric 


ions 


electric  field  region  and  take  advantage  of  the  pressure  gradient  set 


up  by  the  slower  moving  C ° ions.  Since  the  accelerating  electric 


field  is  different  for  the  two  species  Eq.  2 is  not  quite  valid  and  more 


H ions  are  found  at  higher  EA/Z^  than  C ® ions  in  Fig.  4 


Typical  ion  acceleration  time  scale  for  an  initial  electron 


temperature  of  ~ 30  keV  and  scale  length  of  a few  microns  is  a few 


tens  of  picoseconds.  This  relatively  high  electron  temperature  has 


then  to  be  maintained  for  only  a short  time  and  is  therefore  a peak 


temperature  consistent  with  ocher  simulations 


The  multipeak 


structure  on  the  ion  energy  distributions  is  not  reproduced  by  our 


model.  However,  temporal  variations  of  Che  pressure  gradient  on  the 


time  scale  mentioned  above  could  give  bursts  of  ions  of  decreasing 
energies. 

In  summary,  the  electron  pressure  gradient  model  can  account  for 
the  energy  and  the  relative  behavior  of  high  energy  ion  species. 
Furthermore,  it  was  also  shown  that  for  some  plasmas  (CD^  for  example) 
a single  ion  species  description  is  quite  adequate  to  simulate  the 
plasma  expansion  although,  for  others,  a multi-species  description  is 
required.  The  determining  criterion  is  whether  all  the  Ions  have  the 
same  A/Z  during  the  expansion.  A CH^  target  can  never  have  the  same 
A/Z.  Caution  should  also  be  used  with  targets  made  from  higher  Z 
material  (glass,  kt,  etc.)  of  unknown  degree  of  ionization  during  the 
expansion.  The  mixture  of  different  ion  species  results  into  a 
preferential  acceleration  of  the  lower  A/Z  ions  by  the  higher  A/Z  ions. 

The  authors  wish  to  thank  S.  E.  Bodner,  H.  R.  Griem,  F.  S,  Felber 
and  I,  B.  Bernstein  for  their  Interesting  comments  on  this  work  and 
D,  G,  Colombant  for  his  valuable  assistance  in  the  numerical 


calculations. 
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Fig.  2 - (a)  High  energy  ion  distribution  from  a CH2 
target  (b)  oscilloscope  trace  from  a biased  charge 
collector  showing  the  portion  of  the  trace  sampled 
by  the  ion  analyzer. 


THERMAL 


Fig.  3 - Ion  densities,  ion  energies  and  energy  densities  versus 
distance  at  t = 3.  6 r . The  dashed  lines  represent  the  initial  den- 
sity profiles.  The  ion  densities  are  normalized  to  the  initial  H"*" 
plateau  density.  The  non-dimensional  units  are:  r = ?/Cg,  Cg  = 
(hTg^/nOg)! where  T^^  is  the  initial  electron  temperature  and 
C the  initial  density  gradient  scale  length.  The  percentage  in 
parenthesis  gives  the  energy  partition. 
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